No gene by stressful life events interaction on individual differences in adults' self-control.

Background: Difficulty with self-control, or the ability to alter impulses and behavior in a goal-directed way, predicts interpersonal conflict, lower socioeconomic attainments, and more adverse health outcomes. Etiological understanding, and intervention for low self-control is, therefore, a public health goal. A prominent developmental theory proposes that individuals with high genetic propensity for low self-control that are also exposed to stressful environments may be most at-risk of low levels of self-control. Here we examine if polygenic measures associated with behaviors marked by low self-control interact with stressful life events in predicting self-control. Methods: Leveraging molecular data from a large population-based Dutch sample (N = 7,090, Mage = 41.2) to test for effects of genetics (i.e., polygenic scores for ADHD and aggression), stressful life events (e.g., traffic accident, violent assault, financial problems), and a gene-by-stress interaction on self-control (measured with the ASEBA Self-Control Scale). Results: Both genetics (ß =.03 -.04, p <.001) and stressful life events (ß = .11 -.14, p <.001) were associated with individual differences in self-control. We find no evidence of a gene-by-stressful life events interaction on individual differences in adults' self-control. Conclusion: Our findings are consistent with the notion that genetic influences and stressful life events exert largely independent effects on adult self-control. However, the small effect sizes of polygenic scores increases the likelihood of null results. Genetically-informed longitudinal research in large samples can further inform the etiology of individual differences in self-control from early childhood into later adulthood and its downstream implications for public health.

Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

Genome-Wide Association Study of Obsessive-Compulsive Symptoms including 33,943 individuals from the general population.

While 1-2% of individuals meet the criteria for a clinical diagnosis of obsessive-compulsive disorder (OCD), many more (~13-38%) experience subclinical obsessive-compulsive symptoms (OCS) during their life. To characterize the genetic underpinnings of OCS and its genetic relationship to OCD, we conducted the largest genome-wide association study (GWAS) meta-analysis of parent- or self-reported OCS to date (N = 33,943 with complete phenotypic and genome-wide data), combining the results from seven large-scale population-based cohorts from Sweden, the Netherlands, England, and Canada (including six twin cohorts and one cohort of unrelated individuals). We found no genome-wide significant associations at the single-nucleotide polymorphism (SNP) or gene-level, but a polygenic risk score (PRS) based on the OCD GWAS previously published by the Psychiatric Genetics Consortium (PGC-OCD) was significantly associated with OCS (Pfixed = 3.06 × 10-5). Also, one curated gene set (Mootha Gluconeogenesis) reached Bonferroni-corrected significance (Ngenes = 28, Beta = 0.79, SE = 0.16, Pbon = 0.008). Expression of genes in this set is high at sites of insulin mediated glucose disposal. Dysregulated insulin signaling in the etiology of OCS has been suggested by a previous study describing a genetic overlap of OCS with insulin signaling-related traits in children and adolescents. We report a SNP heritability of 4.1% (P = 0.0044) in the meta-analyzed GWAS, and heritability estimates based on the twin cohorts of 33-43%. Genetic correlation analysis showed that OCS were most strongly associated with OCD (rG = 0.72, p = 0.0007) among all tested psychiatric disorders (N = 11). Of all 97 tested phenotypes, 24 showed a significant genetic correlation with OCS, and 66 traits showed concordant directions of effect with OCS and OCD. OCS have a significant polygenic contribution and share genetic risk with diagnosed OCD, supporting the hypothesis that OCD represents the extreme end of widely distributed OCS in the population.

Using Polygenic Scores for Circadian Rhythms to Predict Wellbeing, Depressive Symptoms, Chronotype, and Health.

The association between circadian rhythms and diseases has been well established, while the association with mental health is less explored. Given the heritable nature of circadian rhythms, this study aimed to investigate the relationship between genes underlying circadian rhythms and mental health outcomes, as well as a possible gene-environment correlation for circadian rhythms. Polygenic scores (PGSs) represent the genetic predisposition to develop a certain trait or disease. In a sample from the Netherlands Twin Register (N = 14,021), PGSs were calculated for two circadian rhythm measures: morningness and relative amplitude (RA). The PGSs were used to predict mental health outcomes such as subjective happiness, quality of life, and depressive symptoms. In addition, we performed the same prediction analysis in a within-family design in a subset of dizygotic twins. The PGS for morningness significantly predicted morningness (R2 = 1.55%) and depressive symptoms (R2 = 0.22%). The PGS for RA significantly predicted general health (R2 = 0.12%) and depressive symptoms (R2 = 0.20%). Item analysis of the depressive symptoms showed that 4 out of 14 items were significantly associated with the PGSs. Overall, the results showed that people with a genetic predisposition of being a morning person or with a high RA are likely to have fewer depressive symptoms. The four associated depressive symptoms described symptoms related to decision-making, energy, and feeling worthless or inferior, rather than sleep. Based on our findings future research should include a substantial role for circadian rhythms in depression research and should further explore the gene-environment correlation in circadian rhythms.

Using Instrumental Variables to Measure Causation over Time in Cross-Lagged Panel Models.

Cross-lagged panel models (CLPMs) are commonly used to estimate causal influences between two variables with repeated assessments. The lagged effects in a CLPM depend on the time interval between assessments, eventually becoming undetectable at longer intervals. To address this limitation, we incorporate instrumental variables (IVs) into the CLPM with two study waves and two variables. Doing so enables estimation of both the lagged (i.e., "distal") effects and the bidirectional cross-sectional (i.e., "proximal") effects at each wave. The distal effects reflect Granger-causal influences across time, which decay with increasing time intervals. The proximal effects capture causal influences that accrue over time and can help infer causality when the distal effects become undetectable at longer intervals. Significant proximal effects, with a negligible distal effect, would imply that the time interval is too long to estimate a lagged effect at that time interval using the standard CLPM. Through simulations and an empirical application, we demonstrate the impact of time intervals on causal inference in the CLPM and present modeling strategies to detect causal influences regardless of the time interval in a study. Furthermore, to motivate empirical applications of the proposed model, we highlight the utility and limitations of using genetic variables as IVs in large-scale panel studies.

A multi-omics data analysis workflow packaged as a FAIR Digital Object.

BACKGROUND: Applying good data management and FAIR (Findable, Accessible, Interoperable, and Reusable) data principles in research projects can help disentangle knowledge discovery, study result reproducibility, and data reuse in future studies. Based on the concepts of the original FAIR principles for research data, FAIR principles for research software were recently proposed. FAIR Digital Objects enable discovery and reuse of Research Objects, including computational workflows for both humans and machines. Practical examples can help promote the adoption of FAIR practices for computational workflows in the research community. We developed a multi-omics data analysis workflow implementing FAIR practices to share it as a FAIR Digital Object. FINDINGS: We conducted a case study investigating shared patterns between multi-omics data and childhood externalizing behavior. The analysis workflow was implemented as a modular pipeline in the workflow manager Nextflow, including containers with software dependencies. We adhered to software development practices like version control, documentation, and licensing. Finally, the workflow was described with rich semantic metadata, packaged as a Research Object Crate, and shared via WorkflowHub. CONCLUSIONS: Along with the packaged multi-omics data analysis workflow, we share our experiences adopting various FAIR practices and creating a FAIR Digital Object. We hope our experiences can help other researchers who develop omics data analysis workflows to turn FAIR principles into practice.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Integrative multi-omics analysis of genomic, epigenomic, and metabolomics data leads to new insights for Attention-Deficit/Hyperactivity Disorder.

The evolving field of multi-omics combines data and provides methods for simultaneous analysis across several omics levels. Here, we integrated genomics (transmitted and non-transmitted polygenic scores [PGSs]), epigenomics, and metabolomics data in a multi-omics framework to identify biomarkers for Attention-Deficit/Hyperactivity Disorder (ADHD) and investigated the connections among the three omics levels. We first trained single- and next multi-omics models to differentiate between cases and controls in 596 twins (cases = 14.8%) from the Netherlands Twin Register (NTR) demonstrating reasonable in-sample prediction through cross-validation. The multi-omics model selected 30 PGSs, 143 CpGs, and 90 metabolites. We confirmed previous associations of ADHD with glucocorticoid exposure and the transmembrane protein family TMEM, show that the DNA methylation of the MAD1L1 gene associated with ADHD has a relation with parental smoking behavior, and present novel findings including associations between indirect genetic effects and CpGs of the STAP2 gene. However, out-of-sample prediction in NTR participants (N = 258, cases = 14.3%) and in a clinical sample (N = 145, cases = 51%) did not perform well (range misclassification was [0.40, 0.57]). The results highlighted connections between omics levels, with the strongest connections between non-transmitted PGSs, CpGs, and amino acid levels and show that multi-omics designs considering interrelated omics levels can help unravel the complex biology underlying ADHD.

Longitudinal multi-omics study reveals common etiology underlying association between plasma proteome and BMI trajectories in adolescent and young adult twins.

BACKGROUND: The influence of genetics and environment on the association of the plasma proteome with body mass index (BMI) and changes in BMI remains underexplored, and the links to other omics in these associations remain to be investigated. We characterized protein-BMI trajectory associations in adolescents and adults and how these connect to other omics layers. METHODS: Our study included two cohorts of longitudinally followed twins: FinnTwin12 (N = 651) and the Netherlands Twin Register (NTR) (N = 665). Follow-up comprised 4 BMI measurements over approximately 6 (NTR: 23-27 years old) to 10 years (FinnTwin12: 12-22 years old), with omics data collected at the last BMI measurement. BMI changes were calculated in latent growth curve models. Mixed-effects models were used to quantify the associations between the abundance of 439 plasma proteins with BMI at blood sampling and changes in BMI. In FinnTwin12, the sources of genetic and environmental variation underlying the protein abundances were quantified by twin models, as were the associations of proteins with BMI and BMI changes. In NTR, we investigated the association of gene expression of genes encoding proteins identified in FinnTwin12 with BMI and changes in BMI. We linked identified proteins and their coding genes to plasma metabolites and polygenic risk scores (PRS) applying mixed-effects models and correlation networks. RESULTS: We identified 66 and 14 proteins associated with BMI at blood sampling and changes in BMI, respectively. The average heritability of these proteins was 35%. Of the 66 BMI-protein associations, 43 and 12 showed genetic and environmental correlations, respectively, including 8 proteins showing both. Similarly, we observed 7 and 3 genetic and environmental correlations between changes in BMI and protein abundance, respectively. S100A8 gene expression was associated with BMI at blood sampling, and the PRG4 and CFI genes were associated with BMI changes. Proteins showed strong connections with metabolites and PRSs, but we observed no multi-omics connections among gene expression and other omics layers. CONCLUSIONS: Associations between the proteome and BMI trajectories are characterized by shared genetic, environmental, and metabolic etiologies. We observed few gene-protein pairs associated with BMI or changes in BMI at the proteome and transcriptome levels.

Direct and Indirect Genetic Effects on Aggression.

Background: Family members resemble each other in their propensity for aggression. In twin studies, approximately 50% of the variance in aggression can be explained by genetic influences. However, if there are genotype-environment correlation mechanisms, such as environmental manifestations of parental and sibling genotypes, genetic influences may partly reflect environmental influences. In this study, we investigated the importance of indirect polygenic score (PGS) effects on aggression. Methods: We modeled the effect of PGSs based on 3 genome-wide association studies: early-life aggression, educational attainment, and attention-deficit/hyperactivity disorder (ADHD). The associations with aggression were tested in a within- and between-family design (37,796 measures from 7740 individuals, ages 3-86 years [mean = 14.20 years, SE = 12.03], from 3107 families, 55% female) and in a transmitted/nontransmitted PGS design (42,649 measures from 6653 individuals, ages 3-61 years [mean = 11.81 years, SE = 8.68], from 3024 families, 55% female). All participants are enrolled in the Netherlands Twin Register. Results: We found no evidence for contributions of indirect PGS effects on aggression in either a within- and between-family design or a transmitted/nontransmitted PGS design. Results indicate significant direct effects on aggression for the PGSs based on early-life aggression, educational attainment, and ADHD, although explained variance was low (within- and between-family: early-life aggression R2 = 0.3%, early-life ADHD R2 = 0.6%, educational attainment R2 = 0.7%; transmitted/nontransmitted PGSs: early-life aggression R2 = 0.2%, early-life ADHD R2 = 0.9%, educational attainment R2 = 0.5%). Conclusions: PGSs included in the current study had a direct (but no indirect) effect on aggression, consistent with results of previous twin and family studies. Further research involving other PGSs for aggression and related phenotypes is needed to determine whether this conclusion generalizes to overall genetic influences on aggression.

Longitudinal multi-omics study reveals common etiology underlying association between plasma proteome and BMI trajectories in adolescent and young adult twins.

Background: The influence of genetics and environment on the association of the plasma proteome with body mass index (BMI) and changes in BMI remain underexplored, and the links to other omics in these associations remain to be investigated. We characterized protein-BMI trajectory associations in adolescents and adults and how these connect to other omics layers. Methods: Our study included two cohorts of longitudinally followed twins: FinnTwin12 (N=651) and the Netherlands Twin Register (NTR) (N=665). Follow-up comprised four BMI measurements over approximately 6 (NTR: 23-27 years old) to 10 years (FinnTwin12: 12-22 years old), with omics data collected at the last BMI measurement. BMI changes were calculated using latent growth curve models. Mixed-effects models were used to quantify the associations between the abundance of 439 plasma proteins with BMI at blood sampling and changes in BMI. The sources of genetic and environmental variation underlying the protein abundances were quantified using twin models, as were the associations of proteins with BMI and BMI changes. In NTR, we investigated the association of gene expression of genes encoding proteins identified in FinnTwin12 with BMI and changes in BMI. We linked identified proteins and their coding genes to plasma metabolites and polygenic risk scores (PRS) using mixed-effect models and correlation networks. Results: We identified 66 and 14 proteins associated with BMI at blood sampling and changes in BMI, respectively. The average heritability of these proteins was 35%. Of the 66 BMI-protein associations, 43 and 12 showed genetic and environmental correlations, respectively, including 8 proteins showing both. Similarly, we observed 6 and 4 genetic and environmental correlations between changes in BMI and protein abundance, respectively. S100A8 gene expression was associated with BMI at blood sampling, and the PRG4 and CFI genes were associated with BMI changes. Proteins showed strong connections with many metabolites and PRSs, but we observed no multi-omics connections among gene expression and other omics layers. Conclusions: Associations between the proteome and BMI trajectories are characterized by shared genetic, environmental, and metabolic etiologies. We observed few gene-protein pairs associated with BMI or changes in BMI at the proteome and transcriptome levels.

Environment-by-PGS Interaction in the Classical Twin Design: An Application to Childhood Anxiety and Negative Affect.

One type of genotype-environment interaction occurs when genetic effects on a phenotype are moderated by an environment; or when environmental effects on a phenotype are moderated by genes. Here we outline these types of genotype-environment interaction models, and propose a test of genotype-environment interaction based on the classical twin design, which includes observed genetic variables (polygenic scores: PGSs) that account for part of the genetic variance of the phenotype. We introduce environment-by-PGS interaction and the results of a simulation study to address statistical power and parameter recovery. Next, we apply the model to empirical data on anxiety and negative affect in children. The power to detect environment-by-PGS interaction depends on the heritability of the phenotype, and the strength of the PGS. The simulation results indicate that under realistic conditions of sample size, heritability and strength of the interaction, the environment-by-PGS model is a viable approach to detect genotype-environment interaction. In 7-year-old children, we defined two PGS based on the largest genetic association studies for 2 traits that are genetically correlated to childhood anxiety and negative affect, namely major depression (MDD) and intelligence (IQ). We find that common environmental influences on negative affect are amplified for children with a lower IQ-PGS.

Multimodal Data Integration Advances Longitudinal Prediction of the Naturalistic Course of Depression and Reveals a Multimodal Signature of Remission During 2-Year Follow-up.

BACKGROUND: The ability to predict the disease course of individuals with major depressive disorder (MDD) is essential for optimal treatment planning. Here, we used a data-driven machine learning approach to assess the predictive value of different sets of biological data (whole-blood proteomics, lipid metabolomics, transcriptomics, genetics), both separately and added to clinical baseline variables, for the longitudinal prediction of 2-year remission status in MDD at the individual-subject level. METHODS: Prediction models were trained and cross-validated in a sample of 643 patients with current MDD (2-year remission n = 325) and subsequently tested for performance in 161 individuals with MDD (2-year remission n = 82). RESULTS: Proteomics data showed the best unimodal data predictions (area under the receiver operating characteristic curve = 0.68). Adding proteomic to clinical data at baseline significantly improved 2-year MDD remission predictions (area under the receiver operating characteristic curve = 0.63 vs. 0.78, p = .013), while the addition of other omics data to clinical data did not yield significantly improved model performance. Feature importance and enrichment analysis revealed that proteomic analytes were involved in inflammatory response and lipid metabolism, with fibrinogen levels showing the highest variable importance, followed by symptom severity. Machine learning models outperformed psychiatrists' ability to predict 2-year remission status (balanced accuracy = 71% vs. 55%). CONCLUSIONS: This study showed the added predictive value of combining proteomic data, but not other omics data, with clinical data for the prediction of 2-year remission status in MDD. Our results reveal a novel multimodal signature of 2-year MDD remission status that shows clinical potential for individual MDD disease course predictions from baseline measurements.

Integrative Multi-omics Analysis of Childhood Aggressive Behavior.

This study introduces and illustrates the potential of an integrated multi-omics approach in investigating the underlying biology of complex traits such as childhood aggressive behavior. In 645 twins (cases = 42%), we trained single- and integrative multi-omics models to identify biomarkers for subclinical aggression and investigated the connections among these biomarkers. Our data comprised transmitted and two non-transmitted polygenic scores (PGSs) for 15 traits, 78,772 CpGs, and 90 metabolites. The single-omics models selected 31 PGSs, 1614 CpGs, and 90 metabolites, and the multi-omics model comprised 44 PGSs, 746 CpGs, and 90 metabolites. The predictive accuracy for these models in the test (N = 277, cases = 42%) and independent clinical data (N = 142, cases = 45%) ranged from 43 to 57%. We observed strong connections between DNA methylation, amino acids, and parental non-transmitted PGSs for ADHD, Autism Spectrum Disorder, intelligence, smoking initiation, and self-reported health. Aggression-related omics traits link to known and novel risk factors, including inflammation, carcinogens, and smoking.

Meta-analysis of genome-wide association studies of hoarding symptoms in 27,537 individuals.

Hoarding Disorder (HD) is a mental disorder characterized by persistent difficulties discarding or parting with possessions, often resulting in cluttered living spaces, distress, and impairment. Its etiology is largely unknown, but twin studies suggest that it is moderately heritable. In this study, we pooled phenotypic and genomic data from seven international cohorts (N = 27,537 individuals) and conducted a genome wide association study (GWAS) meta-analysis of parent- or self-reported hoarding symptoms (HS). We followed up the results with gene-based and gene-set analyses, as well as leave-one-out HS polygenic risk score (PRS) analyses. To examine a possible genetic association between hoarding symptoms and other phenotypes we conducted cross-trait PRS analyses. Though we did not report any genome-wide significant SNPs, we report heritability estimates for the twin-cohorts between 26-48%, and a SNP-heritability of 11% for an unrelated sub-cohort. Cross-trait PRS analyses showed that the genetic risk for schizophrenia and autism spectrum disorder were significantly associated with hoarding symptoms. We also found suggestive evidence for an association with educational attainment. There were no significant associations with other phenotypes previously linked to HD, such as obsessive-compulsive disorder, depression, anxiety, or attention-deficit hyperactivity disorder. To conclude, we found that HS are heritable, confirming and extending previous twin studies but we had limited power to detect any genome-wide significant loci. Much larger samples will be needed to further extend these findings and reach a "gene discovery zone". To move the field forward, future research should not only include genetic analyses of quantitative hoarding traits in larger samples, but also in samples of individuals meeting strict diagnostic criteria for HD, and more ethnically diverse samples.

A genetically informed prediction model for suicidal and aggressive behaviour in teens.

Suicidal and aggressive behaviours cause significant personal and societal burden. As risk factors associated with these behaviours frequently overlap, combined approaches in predicting the behaviours may be useful in identifying those at risk for either. The current study aimed to create a model that predicted if individuals will exhibit suicidal behaviour, aggressive behaviour, both, or neither in late adolescence. A sample of 5,974 twins from the Child and Adolescent Twin Study in Sweden (CATSS) was broken down into a training (80%), tune (10%) and test (10%) set. The Netherlands Twin Register (NTR; N = 2702) was used for external validation. Our longitudinal data featured genetic, environmental, and psychosocial predictors derived from parental and self-report data. A stacked ensemble model was created which contained a gradient boosted machine, random forest, elastic net, and neural network. Model performance was transferable between CATSS and NTR (macro area under the receiver operating characteristic curve (AUC) [95% CI] AUCCATSS(test set) = 0.709 (0.671-0.747); AUCNTR = 0.685 (0.656-0.715), suggesting model generalisability across Northern Europe. The notable exception is suicidal behaviours in the NTR, which was no better than chance. The 25 highest scoring variable importance scores for the gradient boosted machines and random forest models included self-reported psychiatric symptoms in mid-adolescence, sex, and polygenic scores for psychiatric traits. The model's performance is comparable to current prediction models that use clinical interviews and is not yet suitable for clinical use. Moreover, genetic variables may have a role to play in predictive models of adolescent psychopathology.

Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention.

Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.

Heritability of Urinary Amines, Organic Acids, and Steroid Hormones in Children.

Variation in metabolite levels reflects individual differences in genetic and environmental factors. Here, we investigated the role of these factors in urinary metabolomics data in children. We examined the effects of sex and age on 86 metabolites, as measured on three metabolomics platforms that target amines, organic acids, and steroid hormones. Next, we estimated their heritability in a twin cohort of 1300 twins (age range: 5.7-12.9 years). We observed associations between age and 50 metabolites and between sex and 21 metabolites. The monozygotic (MZ) and dizygotic (DZ) correlations for the urinary metabolites indicated a role for non-additive genetic factors for 50 amines, 13 organic acids, and 6 steroids. The average broad-sense heritability for these amines, organic acids, and steroids was 0.49 (range: 0.25-0.64), 0.50 (range: 0.33-0.62), and 0.64 (range: 0.43-0.81), respectively. For 6 amines, 7 organic acids, and 4 steroids the twin correlations indicated a role for shared environmental factors and the average narrow-sense heritability was 0.50 (range: 0.37-0.68), 0.50 (range; 0.23-0.61), and 0.47 (range: 0.32-0.70) for these amines, organic acids, and steroids. We conclude that urinary metabolites in children have substantial heritability, with similar estimates for amines and organic acids, and higher estimates for steroid hormones.

DNA methylation in peripheral tissues and left-handedness.

Handedness has low heritability and epigenetic mechanisms have been proposed as an etiological mechanism. To examine this hypothesis, we performed an epigenome-wide association study of left-handedness. In a meta-analysis of 3914 adults of whole-blood DNA methylation, we observed that CpG sites located in proximity of handedness-associated genetic variants were more strongly associated with left-handedness than other CpG sites (P = 0.04), but did not identify any differentially methylated positions. In longitudinal analyses of DNA methylation in peripheral blood and buccal cells from children (N = 1737), we observed moderately stable associations across age (correlation range [0.355-0.578]), but inconsistent across tissues (correlation range [- 0.384 to 0.318]). We conclude that DNA methylation in peripheral tissues captures little of the variance in handedness. Future investigations should consider other more targeted sources of tissue, such as the brain.